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►  Thermal  efficiency  improves  for  certain  range  of  pin  geometry. 

►  Thermal  stress  in  thermoelectric  module  varies  for  different  configuration  of  pin  geometry. 

►  Maximum  stress  occurs  around  the  edges  of  the  pins  in  the  region  close  to  the  hot  plate. 

►  The  pin  geometry  with  RA  =  2  results  in  slightly  low  thermal  stress  in  the  pin. 
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Thermodynamics  and  thermal  stress  analysis  of  thermo-electric  power  generator  is  carried  out.  The 
influence  of  device  geometry  on  thermal  stress,  thermal  efficiency  and  output  power  is  examined.  The 
finite  element  method  is  incorporated  to  predict  temperature  and  stress  fields  in  the  thermo  electric 
device.  It  is  found  that  thermal  efficiency  improves  for  certain  geometric  configuration  of  the  device.  In 
this  case;  the  maximum  thermal  stress  developed  in  the  pin  reduces  slightly  indicating  improved  life 
expectation  of  the  device. 

©  2012  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Thermoelectric  power  generator  is  one  of  the  solid-state  devices 
converting  thermal  energy  into  an  electrical  energy.  Thermoelectric 
generators  operate  in  between  the  high  and  low  temperature 
sources  and  the  efficiency  of  the  device  increases  with  increasing 
temperature  difference  between  the  sources.  However,  thermal 
stress  developed  within  the  device  limits  temperature  difference  in 
the  practical  applications  due  to  the  shorting  of  the  life  cycle  of  the 
device.  Moreover,  material  failure  due  to  high  stress  induced 
cracking  prevents  further  operations  of  the  device  with  expected 
performance.  Consequently,  investigation  into  thermal  stress 
development  in  the  thermoelectric  device  becomes  essential. 

The  efficiency  of  thermoelectric  power  generator  under  thermal 
cyclic  loading  was  investigated  by  Hatzikraniotis  et  al.  [1].  They 
indicated  that  cyclic  thermal  loading  of  the  device  resulted  in 
reduction  in  Seebeck  current  and  increased  in  electric  and  thermal 
resistivity.  Thermal  expansion  of  high  temperature  thermoelectric 
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materials  was  investigated  by  the  Ravi  et  al.  [2].  They  showed  that 
thermal  expansion  coefficient  of  the  thermoelectric  materials  was 
the  critical  parameter  influencing  stress  levels  in  the  device.  A  new 
method  for  thermoelectric  power  generation  using  large  area  p-n- 
junction  was  presented  by  Span  et  al.  [3].  Their  findings  revealed 
that  the  new  design  of  the  device  eliminated  the  thermal  stress 
between  active  elements,  such  as  pins,  and  joining  material. 
Thermo-mechanical  modeling  of  large  area  thermoelectric  gener¬ 
ators  was  carried  out  by  Turenne  et  al.  [4].  They  indicated  that  the 
maximum  stress  was  increased  in  the  thermoelectric  device  close 
to  the  edge  of  the  pins.  Thermo-mechanical  analysis  of  thermo¬ 
electric  power  generator  was  carried  out  by  Soto  and  Venkatasu- 
bramanian  [5].  They  showed  that  due  to  the  large  temperature 
gradients  as  well  as  mismatch  in  thermal  expansion  coefficients, 
high  stress  levels  were  formed  in  the  thermoelectric  device. 
Thermal  stress  analysis  for  thermoelectric  generator  was  carried 
out  by  Nakatani  et  al.  [6].  The  findings  revealed  aluminum  and 
molybdenum  coatings  formed  by  atmospheric  plasma  spraying 
reduce  the  thermal  stain  at  the  interface  of  the  pins.  Thermal 
expansion  in  the  thermoelectric  generators  was  examined  by 
Hilkage  at  [7].  They  showed  that  the  grain  size  in  the  active 
elements  influenced  significantly  thermal  expansion  coefficient  of 


684 


AS.  Al-Merbati  et  al.  /  Applied  Thermal  Engineering  50  (2013)  683-692 


the  device  pins.  Thermal  efficiency  of  the  topping  cycle  with  and 
without  presence  of  thermoelectric  generator  was  investigated  by 
Sahin  et  al.  [8].  They  demonstrated  that  thermal  efficiency  of  the 
topping  cycle  becomes  slightly  higher  than  that  of  the  same  system 
without  the  presence  of  the  thermoelectric  generators  for  a  certain 
combination  of  operating  and  thermoelectric  device  parameters. 
Numerical  modeling  and  design  of  thermoelectric  cooling  systems 
were  carried  out  by  Jou  [9].  He  introduced  three  approaches  to 
obtain  the  internal  parameters  of  a  given  thermoelectric  module  for 
designing  of  a  thermoelectric  cooling  system.  Thermoelectric 
cooling  of  microelectronic  circuits  and  waste  heat  recovery  were 
examined  by  Gould  et  al.  [10].  They  demonstrated  that  maximum 
electrical  power  could  be  generated  by  the  thermoelectric  module 
when  a  water  cooled  heat  sink  was  used  as  the  secondary  heat  sink, 
as  this  produced  the  greatest  temperature  difference  between  both 
sides  of  the  module.  The  irreversible  analysis  of  a  multi-element 
thermoelectric  generator  system  was  carried  out  by  Xiao  et  al. 
[11  ].  They  indicated  found  that  the  existence  of  the  irreversible  heat 
convection  process  caused  a  large  loss  of  heat  exergy  in  the  ther¬ 
moelectric  genetor  system  and  using  thermoelectric  generators  for 
low-grade  waste  heat  recovery  had  a  promising  potential. 

Although  considerable  research  studies  were  carried  out  to 
examine  the  thermodynamic  performance  of  the  thermoelectric 
device,  thermal  stress  developed  due  to  temperature  gradients  is  not 
given  enough  attention.  In  addition,  the  location  of  the  maximum 
stress  level  in  the  thermoelectric  device  and  stress  the  stres  intesnity 
minimization  through  altering  geometric  configuration  of  device  pins 
are  left  obscure  in  the  previous  studies.  Consequently,  in  the  present 
study,  thermal  efficiency  and  thermal  stress  developed  in  the  ther¬ 
moelectric  generator  are  investigated.  The  temperature  and  stress 
fields  in  the  device  during  the  operation  are  stimulated  using  the  finite 
element  method.  The  influence  of  device  pin  geometry  on  the  stress 
levels  is  examined  and  improvement  in  thermodynamic  efficiency  due 
to  geometric  configuration  of  the  thermoelectric  device  is  presented. 


2.  Mathematical  analysis 

The  analysis  pertinent  to  thermoelectric  generator  is  divided 
into  two-sub  sections  including  the  thermodynamic  analysis  and 
thermal  stress  formulations. 

2.1.  Thermodynamics  analysis 


The  efficiency  of  the  trapezoidal  pin  thermoelectric  power 
generator,  which  is  shown  in  Fig.  1  is  given  as  [12]: 


_ ERl _ 

alTH+K(TH-TL)  -1 12R 


(1) 


where  I<  is  the  thermal  conductance  and  R  is  the  electrical  resis¬ 
tivity  of  the  thermoelectric  generator. 

The  current  /  is  a  function  of  Seebeck  coefficient  a  =  ap  -  an, 
the  upper  and  lower  junction  (plate)  temperatures  (TH  and  TL),  the 
electrical  resistance  R  and  the  external  load  resistance  RL  as 


/  = 


«(Th  -  T 


V  = 


Rl  +  R  ’  (2) 

Substituting  equation  (2)  in  equation  (1)  the  efficiency  becomes 

a2(TH  -  Tl)Rl 


K(Rl  +  R)2+a2TH(RL  +  R)  -  L2(TH  -  TL)R 


(3) 


The  cross-sectional  area  of  the  trapezoidal  leg  of  the  thermo¬ 
electric  generator  shown  in  Fig.  2  is 


Fig.  1.  Schematic  view  of  a  thermoelectric  power  generator  and  pin  configurations. 


A(x)  =  —  —  —x  +  Al 


(4) 


where  Al  is  the  cross  sectional  area  of  the  bottom  side  of  the  leg  and 
Ah  is  that  of  the  top  side.  L  is  the  height  of  the  leg.  Defining  an  area 
ratio  Ra  =  Ah/Al,  the  cross-sectional  area  of  the  leg  can  be  written  as 


A{x) 


1+2 


Ra  —  1 
Ra+ 1 


(5) 


where  Ao  is  the  cross-sectional  area  of  the  uniform  leg. 

The  heat  transfer  rate  through  the  leg  along  x  is  given  by: 

Q=-kA(x)^  (6) 

After  assuming  a  steady  heating  situation  and  isolated  leg 
surfaces,  equation  (6)  can  be  re-arranged  as 
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+§r+" 


dr 


(7) 


A  /ft a 


- 1 


I<  =  2(/<p  +  /<n) 


L  V^A  +  1 

ln(ftA) 


(10) 


Making  use  of  equation  (5)  in  equation  (7),  and  performing  the 
integration 


a  - 


2k—  |^A~  1 
t  l«A+1 

ln(RA) 


(I'm  -  Th 


(8) 


where  kp  and  kn  are  the  thermal  conductivities  of  the  p-type  and  n- 
type  legs,  respectively. 

On  the  other  hand,  the  overall  electrical  resistance  of  the  leg  can 
be  written  as 


R\eg  — 


L. 

/ 


dx 

kA(x) 


(11) 


Equation  (8)  indicates  that  the  overall  thermal  conductance  of 
the  trapezoidal  leg  is 


Substituting  A(x)  from  equation  (5)  and  performing  the  inte¬ 
gration,  the  overall  electrical  resistance  is  obtained  as 


2k 


A  /ft A 


- 1 


Kleg  — 


L  VftA  +  1 

ln(ftA) 


(9) 


ftleg  — 


2  k( 


A  /ft A 


- 1 


L  VftA  +  1 


ln(ftA) 


(12) 


Considering  the  two  legs  in  Fig.  1,  the  total  thermal  conductance 
of  the  thermoelectric  generator  can  be  written  as 


Similarly,  considering  the  two  legs  the  total  electrical  resistance 
of  the  thermoelectric  generator  becomes 


3-Dimensional  view  of  thermoelectric  generator 
Maximum  Stress  Line  Centerline 
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I  i 
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Fig.  3.  Three-dimensional  view  of  thermoelectric  generator  and  its  dimensions  used  in  the  simulations. 
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R  = 


RA_l\ln(RA) 

RTT v 


ke.p  +  k, 


e,n 


2k  k  /Ra  1 

^<'e,p Ke,n~i~  <  - 


\Ra 


In (Ra) 


(13) 


where  kep  and  ken  are  the  electrical  conductivities  of  the  p-type  and 
n-type  legs,  respectively.  Substituting  equations  (10)  and  (13)  in 
equation  (3)  the  efficiency  of  the  thermoelectric  generator  can  be 
written  in  dimensionless  form  as 


£  =  (rk  +  l)  [KA  -  1 

K0  in (Ra)  [Ra  +  1. 

and 


R 
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In  (Ra) 

Ro  ~ 
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RL  1  K  „ 

R0  +  2R0(l 


■0) 


(14) 


where 


Table  2 

Properties  of  Bi2Te3  [15-17]. 


T 

6  =  (Temperature  ratio) 

(15) 

k 

rk  =  (Thermal  conductivity  ratio) 

kn 

(16) 

rke  m  (Electrical  conductivity  ratio) 

ke.n 

(17) 

(The  figure  of  merit  based  on  the  average  temperature)  (18) 
I<0  =  (Reference  thermal  conductance)  (19) 


Temperature  (I<) 

Thermal  conductivity  (W/m-K) 

325 

0.93 

375 

0.9 

425 

0.91 

475 

0.95 

525 

1.1 

Specific  heat  (J/K-kg) 

Density  (kg/m~3) 

Yield  stress  (Pa) 

154.4 

7740 

1.12E+08 

Thermal  expansion  (/ 1<) 

Temperature  (I<) 

8.00E-06 

297 

1.01E-05 

304.3 

1.21E-05 

365 

1.24E-05 

451 

1.32E-05 

613 

1.33E-05 

793 

1.41E-05 

864 

Poisson’s  ratio 

Young’s  modulus  (Pa) 

Temperature  (I<) 

0.23 

6.5E  +  010 

200 

0.23 

6.3E  +  010 

300 

0.23 

6.2E  +  010 

400 

0.23 

6.0E  +  010 

500 

0.23 

5.9E  +  010 

600 

and 

Rq  =  ~z  — (Reference  electrical  resistivity).  (20) 

AoKe^n 

Thus,  the  overall  thermal  conductance  and  overall  electrical 
resistivity  can  be  written  in  dimensionless  form,  respectively,  as 


Table  1 

Properties  of  Bi3Te2  used  for  the  efficiency 
calculation. 
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Fig.  4.  Validation  of  present  simulations  with  the  previous  study  [11]. 
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Fig.  5.  Thermal  efficiency  variation  with  RA  (pin  area  ratio  of  at  hot  junction  to  cold 
junction,  RA  =  AHjAL. 


2.2.  Thermal  stress  formulations 

Fig.  3  shows  the  schematic  view  of  three-dimensional  thermo¬ 
electric  generator.  The  thermoelectric  module  has  inhomogeneous 
thermal  and  mechanical  properties  in  the  z-direction  since  the 
properties  used  in  this  study  are  temperature  dependent.  The 
transient  heat  conduction  equation  considered  is: 


0 

Ox 


7  /  \ 

klz>lTx 


or 

ay. 


or 

=  CpP  ot 


(23) 


The  coupled  thermal  stress  analysis  require  to  identify  the 
displacement— strain  relations,  which  are  expressed  in  dimen¬ 
sionless  form  as  follows  [13]: 


0  u  _  dv  0  w 

=  W  £yy  =  W  £ZZ  =  9Z 


(24) 
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Fig.  6.  (a)Temperature  distribution  along  the  centerline  and  the  maximum  stress  line 
for  RA  =  0.5.(b)Temperature  gradient  along  the  centerline  and  the  maximum  stress  line 
(25)  for  Ra  =  0.5. 


An  exact  implementation  of  Newton’s  method  involves 
a  nonsymmetrical  Jacobian  matrix  which  is  stress-strain  relation 
in  dimensionless  form  as  is  illustrated  in  the  following  matrix 
representation  of  the  coupled  equations  [13]: 


3.  Numerical  simulations 

The  thermal  efficiency,  which  is  presented  in  equation  (14)  is 
simulated  using  the  data  in  Table  1.  For  the  thermal  stress  analysis 
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(26) 


Solving  this  system  of  equations  requires  the  use  of  the 
unsymmetrical  matrix  storage  and  solution  scheme.  Furthermore, 
the  mechanical  and  thermal  equations  are  solved  simultaneously. 


ABAQUS  finite  element  code  [14]  is  used  to  simulate  temperature 
and  stress  fields  in  the  thermoelectric  generator.  Table  2  gives  the 
data  used  in  the  simulations. 
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In  the  thermal  stress  analysis  simulations,  it  is  considered  that 
the  thermoelectric  generator  consists  of  ceramic  substrate,  copper 
plate,  tin-Lead  solder,  and  thermoelectric  pins  (Fig.  3).  It  is  also 
considered  that  the  thickness  of  copper  plate  is  0.12  mm,  the 
thickness  of  solder  is  0.04  mm,  and  the  thickness  of  ceramic 
substrate  is  0.34  mm.  The  size  of  the  thermoelectric  generator  pins 
is  3  mm  x  3  mm  x  3  mm  (Fig.  3).  In  the  thermal  stress  simulations, 
it  is  assumed  that  the  thermoelectric  pins  are  made  from  Bi2Te3  and 
that  there  is  no  difference  in  properties  as  a  function  of  position,  so 
that  the  thermoelectric  and  mechanical  properties  are  function  of 
temperature  only.  The  thermal  conductivity  k  (T),  coefficient  of 
liner  thermal  expansion  a  (T),  specific  heat  capacity  cp  (T),  and 
modulus  of  elasticity  E  (T)  are  the  function  of  temperature. 

To  validate  the  predictions  of  thermal  stress  analysis,  simulations 
are  extended  to  include  the  case,  which  was  presented  in  the  early 
study  [15].  In  this  case,  simulations  conditions  are  re-set  to  be  in  line 


with  the  previous  study  [15]  the  data  used  in  the  previous  study  are 
incorporated  in  the  simulations.  The  thermal  stress  predictions  of 
the  present  study  and  the  results  of  the  previous  study  is  shown  in 
Fig.  4.  It  is  evident  that  both  results  are  in  good  agreement.  Therefore, 
the  model  used  in  the  present  study  is  justifiable  to  proceed  with  the 
predictions  of  thermal  stress  fields  in  the  thermoelectric  generator. 

4.  Results  and  discussion 

Thermodynamic  analysis  and  thermal  stress  field  developed  in 
the  thermoelectric  generator  are  investigated  and  the  influence  of 
pin  geometric  configuration  on  thermal  efficiency  and  thermal 
stress  levels  are  examined. 

Fig.  5  shows  efficiency  variation  with  the  geometric  parameter 
Ra,  which  is  the  area  ratio  of  the  pin  cross-section  at  cold  junction 


a 


514.00 

49S.91 

477.83 

4S9.74 

441.65 

423.57 

405.48 

387.39 

369.31 

351.22 

333.13 

315.04 

296.96 


OD8:  R_K«K.odb  4.10-1  W.d  F*b  13  10:30:37  Ar«b  St*f»4*rd  Tim*  2012 


Step:  heating 

Increment  85:  Step  Time  -  80.00 

Primary  Ver:  NTH 


b 


•►1. 860409 
41.710409 
41.550409 
41.400409 
41.240409 
41.090409 
49.330408 
47.780408 
46.220408 
44.670408 
43.120408 
41.570408 
41.900406 


OOB:  R_h*K.*db  Ab«eus/$t*ftd»r4  4.10*  1  W*d  F«b  15  10:50:37  Ar«b  Tim*  201  2 


Step:  heating 

Increment  85:  Step  Time  -  80.00 

Primary  Var:  S,  Mtses 


Fig.  7.  (a)Three-dimensional  temperature  distribution  in  the  thermoelectric  generator  for  RA  =  0.5.(b)Three-dimensional  thermal  stress  distribution  in  the  thermoelectric  generator 
for  Ra  =  0.5. 
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to  hot  junction  {Ra  =  Ah/Al).  It  can  be  observed  that  thermal  effi¬ 
ciency  of  the  thermoelectric  generator  increases  for  RA  values 
different  than  R&  =  1  and  the  minimum  efficiency  occurs  for  Ra  =  1. 
This  indicates  that  parallel  sided  pins  for  the  thermoelectric 
generator  are  not  always  the  optimum  design  for  the  maximum 
efficiency.  Consequently,  the  pin  design  including  either  RA  <  0.5 
and  Ra  <2  results  in  improved  thermal  efficiency.  Although  the 
improvement  in  thermal  efficiency  is  low,  small  improvement  of 
thermal  efficiency  for  one  device  adds  to  the  improvement  of  the 
overall  thermal  efficiency  of  all  devices  in  series,  which  are  used  for 
the  practical  applications. 

Fig.  6(a)  shows  temperature  distribution  along  the  two  lines  in 
the  y-axis  for  Ra  =  0.5.  The  centerline  represents  the  line  passing 
through  the  pin  center  along  the  y-axis  (Fig.  3)  and  maximum  stress 
line  corresponds  to  the  line  passing  through  the  maximum  stress 
locations  along  the  y-axis  (Fig.  3).  Temperature  distribution  is  pre¬ 
sented  when  the  temperature  in  the  pin  reaches  to  a  steady  state 
distribution.  It  should  be  noted  that  the  transient  heating  takes  place 
only  0.3  s  after  the  heating  is  initiated  and  temperature  field  settles 
through  heat  in  the  thermoelectric  device  including  the  pins.  The 
temperature  decays  gradually  along  the  pin  height  for  both  lines  in 
the  y-axis.  Flowever,  the  temperature  gradients  differ  totally  along 
these  lines,  which  are  evident  from  Fig.  6(b),  in  which  the  temper¬ 
ature  gradient  along  the  y-axis  is  shown.  The  temperature  gradient 
attains  high  values  in  the  region  close  to  the  hot  junction  of  the  pin, 
which  is  more  pronounced  along  the  maximum  stress  line.  Flowever, 
the  temperature  gradient  reduces  toward  the  cold  junction  of  the 
pin.  When  comparing  the  temperature  gradients  along  both  lines;  it 
is  evident  that  the  temperature  gradient  is  higher  in  hot  junction  of 
the  pin  along  the  maximum  stress  line  as  compared  to  that  of  the 
center  line.  Flowever,  opposite  is  true  in  the  region  toward  the  cold 
junction.  Consequently,  large  changes  in  the  temperature  gradient 
results  development  of  high  strain  in  the  pin.  Fig.  7(a)  shows 
temperature  contours  in  the  thermoelectric  generator  and  the  pins. 
It  is  evident  that  temperature  remains  almost  uniform  along  the  x 
and  z-axes.  The  large  variation  in  the  temperature  occurs  along  the 
y-axis  between  the  hot  and  cold  junctions. 

Fig.  8  shows  thermal  stress  variation  along  the  centerline  and 
the  maximum  stress  line  in  the  y-direction  for  Ra  =  0.5.  It  is  evident 
that  thermal  stress  attains  high  values  in  the  region  of  the  hot 
junction  which  is  particularly  true  along  the  maximum  stress  line. 
Thermal  stress  reduces  sharply  in  the  region  close  to  the  hot 
junction  in  the  pin  and  its  decay  become  gradual  toward  the  cold 
junction.  The  sharp  decay  of  thermal  stress  is  attributed  to  attain¬ 
ment  of  high  temperature  gradients  in  this  region.  Fig.  7(b)  shows 
thermal  stress  contours  in  the  thermoelectric  device  and  its  pins. 
High  stress  region  occurs  locally  in  the  pin,  particularly  at  the  edges 
of  the  pin  where  pin  is  attached  to  high  temperature  plate.  The 
attainment  of  high  stress  is  because  of  one  on  all  of  the  following 
reasons:  i)  high  temperature  gradient  developed  in  this  region 
gives  to  high  thermal  stress  levels,  and  ii)  the  difference  in  thermal 
expansion  coefficients  due  to  pin  and  the  hot  plate  which  generates 
high  stress  levels  at  the  interface  location  between  the  hot  plate 
and  the  pin.  Moreover,  low  stress  region  in  the  pin  extends  toward 
the  cold  junction  region. 

Fig.  9(a)  shows  temperature  distribution  along  the  centerline 
and  the  maximum  stress  line  along  the  y-axis  for  Ra  =  1.  It  should 
be  noted  that  Ra  =  1  corresponds  to  rectangular  parallel  pin 
geometry.  Temperature  variation  between  the  hot  and  the  cold 
junctions  along  the  both  lines  is  gradual,  provided  that  the 
temperature  gradient  changes  significantly  in  the  hot  plate  region 
for  both  lines  (Fig.  9b).  In  this  case  the  temperature  gradient  attains 
higher  values  in  the  vicinity  of  the  hot  plate,  where  y  =  0,  for  the 
maximum  stress  line  as  compared  to  that  corresponding  to  the 
centerline.  The  high  temperature  gradient  generates  high  strain  in 
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Fig.  8.  Thermal  stress  distribution  along  the  centerline  and  the  maximum  stress  line 
for  Ra  =  0.5. 

this  region.  Fig.  10(a)  shows  temperature  contours  in  the  thermo¬ 
electric  generator  and  the  pins.  It  is  evident  that  temperature 
variation  along  x  and  z-axes  are  almost  uniform.  Flowever, 
temperature  variation  along  the  y-axis  is  large  due  to  the  presence 
of  hot  and  cold  plates  across  the  pins. 

Fig.  11  shows  thermal  stress  distribution  along  the  centerline 
and  the  maximum  stress  line  in  the  y-axis  for  Ra  =  1.  Thermal  stress 


y  (m) 


x  (m) 

Fig.  9.  (a)Temperature  distribution  along  the  centerline  and  the  maximum  stress  line 
for  RA  =  l.(b)Temperature  gradient  distribution  along  the  centerline  and  the  maximum 
stress  line  for  RA  =  1. 
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varies  significantly  along  the  both  lines  and  stress  reduces  slightly 
in  the  region  next  to  the  hot  plate  and  increases  gradually  with 
increasing  distance  toward  the  cold  plate.  However,  thermal  stress 
reduces  sharply  in  the  region  close  to  the  hot  plate  along  the 
maximum  stress  line  and  the  decay  becomes  gradual  toward  the 
cold  plate.  The  attainment  of  high  stress  along  maximum  stress  line 
is  due  to  thermal  expansion  of  the  pin  and  high  temperature 
gradient  developed  in  the  region  close  to  the  hot  plate.  Fig.  10(b) 
shows  thermal  stress  counters  in  the  thermoelectric  generator  and 
its  pins.  It  can  be  observed  that  high  stress  region  is  locally  devel¬ 
oped  toward  the  pin  edges,  particularly  in  the  region  next  to  the  hot 
plate.  However,  thermal  stress  drops  significantly  toward  to  the 
cold  plate  due  to  attainment  of  the  low  temperature  gradients  in 
this  region. 

Fig.  12(a)  shows  temperature  distribution  along  the  centerline 
and  the  maximum  stress  line  in  the  y-axis  for  Ra  =  2.  Temperature 
decays  gradually  along  the  both  lines  form  the  high  temperature 
plate  to  low  temperature  plate  as  similar  those  shown  in  Figs.  6  and 
9(a).  However,  temperature  along  the  maximum  stress  line  attains 
slightly  higher  values  than  that  of  the  centerline.  This  is  associated 
with  the  heat  conduction  along  the  thermoelectric  pins.  Moreover, 
the  temperature  gradient  increases  gradually  toward  the  cold  plate, 
which  is  associated  with  the  gradual  decay  of  temperature  along 
the  both  lines  (Fig.  12(b)).  The  temperature  gradient  along  the 
maximum  stress  line  is  slightly  lower  than  that  of  along  the 
centerline.  This  indicates  the  attainment  of  relatively  lower  thermal 
strain  along  this  line.  Fig.  13(a)  shows  temperature  contours  in  the 
thermoelectric  device  and  in  the  pins.  Temperature  variation  is 
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Fig.  11.  Thermal  stress  distribution  along  the  centerline  and  the  maximum  stress  line 
for  Ra  =  1. 


almost  uniform  along  the  x  and  z-axes  as  similar  to  those  obtained 
for  the  cases  Ra  =  0.5  and  Ra  =  1. 

Fig.  14  shows  von  Mises  Stress  distribution  along  the  centerline 
and  the  maximum  stress  line  in  the  y-axis  direction  for  Ra  =  2. 
Thermal  stress  behavior  is  similar  to  that  shown  in  Fig.  11  for  Ra  =  1, 
provided  that  the  maximum  stress  at  x  =  0  is  less  than  that  of 
Ra  =  1.  This  is  attributed  to  the  differences  in  the  temperature 
gradient  behavior  particularly  in  the  region  close  to  the  hot  plate. 
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Fig.  10.  (a)Three-dimensional  temperature  distribution  in  the  thermoelectric  gener¬ 
ator  for  Ra  =  l.(b)Three-dimensional  thermal  stress  distribution  in  the  thermoelectric 
generator  for  RA  =  1. 


Fig.  12.  (a)Temperature  distribution  along  the  centerline  and  the  maximum  stress  line 
for  Ra  =  2.(b)Temperature  gradient  distribution  along  the  centerline  and  the  maximum 
stress  line  for  RA  =  2. 
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Fig.  13.  (a)Three-dimensional  temperature  distribution  in  the  thermoelectric  gener¬ 
ator  for  Ra  =  2.(b)Three-dimensional  thermal  stress  distribution  in  the  thermoelectric 
generator  for  RA  =  2. 


Fig.  13(b)  shows  thermal  stress  contours  in  the  thermoelectric 
generator  and  in  the  pins.  Thermal  stress  is  high  in  the  region  close 
to  the  hot  plate  and  reduces  sharply  in  the  region  close  to  the  cold 
plate.  The  presence  of  high  stress  centers  in  the  region  of  the  pin 
edges  is  evident.  This  is  attributed  to  the  attainment  of  high 
temperature  gradient  and  thermal  expansion  coefficient  differ¬ 
ences  between  the  hot  plate  and  the  pin  material.  Table  3  gives  the 
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Fig.  14.  Thermal  stress  distribution  along  the  centerline  and  the  maximum  stress  line 
for  Ra  =  2. 


Table  3 

Comparison  of  the  maximum  and  the  minimum  stresses  predicted  from  the  present 
study  and  obtained  from  the  previous  study  [15]. 


m 

o' 

II 

< 

04 

Ra=1 

Ra  =  2 

Max.  stress  (GPa) 

1.12 

0.969 

0.912 

maximum  values  of  thermal  stress  developed  in  the  thermoelectric 
generator  when  the  device  operation  reaches  to  steady  state 
condition.  The  maximum  thermal  stress  for  Ra  =  2  is  the  minimum 
among  the  other  cases.  This  suggests  that  changing  the  design  of 
the  pins  of  the  thermoelectric  power  generators  form  parallel  to 
trapezoidal  geometry,  the  maximum  thermal  stress  developed  can 
be  reduced.  Although  the  reduction  is  small,  the  life  extension  of 
the  thermoelectric  device  is  possible. 

5.  Conclusion 

Thermodynamic  and  thermal  stress  analysis  of  thermoelectric 
generator  is  carried  out  for  various  device  pin  centrifugations.  In  the 
analysis,  the  properties  of  Bi2Te3  thermoelectric  material  are  used.  It  is 
found  that  thermal  efficiency  of  the  thermoelectric  generator  improves 
for  Ra  <  0.5  and  Ra  <  2.  Although  this  improvement  is  small,  its  effect 
becomes  high  for  the  series  operation  of  many  thermoelectric  devices. 
Temperature  variation  in  the  x,y  plane  does  not  alter  significantly;  in 
which  case,  almost  uniform  temperature  distribution  takes  place  along 
the  pins.  Temperature  variation  along  the  y-axis  results  in  sharp  decay 
of  temperature  in  the  region  close  to  the  hot  plate.  This  gives  rise  to 
attainment  of  high  temperature  gradients  in  this  region.  Thermal  stress 
developed  in  the  thermoelectric  module  and  its  pins  varies  for  different 
configuration  of  pin  geometry  and  the  maximum  stress  occurs  around 
the  edges  of  the  pins,  which  is  more  produced  in  the  region  close  to  the 
hot  plate.  The  attainment  of  high  stress  levels  in  this  region  is  attributed 
to  the  high  temperature  gradients  and  the  mismatch  of  thermal 
expansion  coefficients  of  the  hot  plate  and  the  device  pins.  The  influ¬ 
ence  of  pin  geometry  on  temperature  gradients  and  thermal  stress  is 
notable.  In  this  case  the  pin  with  Ra  =  2  results  in  slightly  low  thermal 
stress  in  the  pin.  This  indicates  that  changing  the  pin  geometry 
improves  temperature  variation  in  the  thermoelectric  generator  while 
suppressing  the  maximum  stress  levels  in  the  pin.  Consequently,  life 
expectation  of  the  device  can  be  improved  with  geometric  design  of  the 
device  pins.  In  addition,  for  pin  geometric  configuration  RA  =  2,  thermal 
efficiency  of  the  device  also  improves. 
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